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ABSTRACT 



Galaxy clusters are being assembled today in the most energetic phase of hierarchical structure formation which manifests itself in powerful 
shocks that contribute to a substantial energy density of cosmic rays (CRs). Hence, clusters are expected to be luminous gamma- ray emitters 
since they also act as energy reservoirs for additional CR sources, such as active galactic nuclei and supernova-driven galactic winds. To detect 
the gamma-ray emission from CR interactions with the ambient cluster gas, we conducted the deepest to date observational campaign targeting 
• a galaxy cluster at very high-energy gamma-rays and observed the Perseus cluster with the MAGIC Cherenkov telescopes for a total of ~ 85 h 
" of effective observing time. This campaign resulted in the detection of the central radio galaxy NGC 1275 at energies E > 100 GeV with a very 
' steep energy spectrum. Here, we restrict our analysis to energies E > 630 GeV and detect no significant gamma-ray excess. This constrains the 
. ^ , average CR-to-thermal pressure ratio to be < 1-2%, depending on assumptions and the model for CR emission. Comparing these gamma-ray upper 
> limits to models inferred from cosmological cluster simulations that include CRs constrains the maximum CR acceleration efficiency at structure 
» I ■ formation shocks to be < 50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into 
' the outer cluster regions. Finally, we derive lower limits on the magnetic field distribution assuming that the Perseus radio mini-halo is generated 
by secondary electrons/positrons that are created in hadronic CR interactions: assuming a spectrum of around TeV energies as implied by 
cluster simulations, we limit the central magnetic field to be > 4-9 fiG, depending on the rate of decline of the magnetic field strength toward larger 
radii. This range is well below field strengths inferred from Faraday rotation measurements in cool cores. Hence, the hadronic model remains a 
plausible explanation of the Perseus radio mini-halo. 

Key words. Gamma rays: galaxies: clusters: individual: Perseus 



1. Introduction 

In the hierarchical model of structure formation, clusters of 
galaxies are presently forming through mergers of smaller 
galaxy groups, and assembling masses up to a few times lO'^'M©. 
During these merger events, enormous amounts of energy — of 
the order of the final gas binding energy fibind ~ 3 x (10^'- 



10 ) erg — are released through collisionless shocks and tur- 
bulence. This energy is dissipated on a dynamical timescale of 
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Tdyn ~ 1 Gyr (see IVoitjl2005L for a review). Hence, the corre- 
sponding energy dissipation rates are L ~ (10'*^-10^^) erg s"^ 
If only a small fraction of this energy were converted into non- 
thermal particle populations, the associated emission should be 
detectable in the gamma-ray regime and can potentially be used 
to decipher the short- and long-term history of structure forma- 
tion. 

Many galaxy clusters show large scale diffuse synchrotron 
radio emission in the form of so-called radio (mini-)halos which 
proves the existence of magnetic fields and relativistic electrons 
perme ating the intra-cluster medium (ICM) (e.g., Kempner et al] 
l2004l:iFeretti et alj2004 ). Through this synchrotron emission we 
can identify the sites of acceleration and injection of the rela- 
tivistic electrons into the ICM: (1) merger shocks propagating 
toward the cluster outskirts, observed as giant radio relics, (2) 
active galactic nuclei (AGN) in clusters, observed through ra- 
dio jets and bubbles, and (3) galactic winds and ram pressure 
stripping are often accompanied by synchrotron emission. In 
analogy with shocks within our Galaxy, such as those in su- 
pernova remnants, galaxy clusters should also be acceleration 
sites for relativistic protons and heavier relativistic nuclei. Due 
to their higher masses compared with the electrons, protons and 
nuclei are accelerated more efficiently to relativistic energies and 
are expected to show a ratio of the spectral energy flux of cos- 
mic ray (CR) protons to CR electrons above 1 GeV of about 
100 (as it is obse rved in our Galaxy between 1-10 GeV, see 
ISchhckeiseil2002l) . CR protons also have radiative cooling times 
that are larger than the corresponding cooling times of CR elec- 
trons by the square of the mass ratio, (mp/mg)^, and hence can 
accumulate for the Hubble time in a galaxy cluster dVoIk et al] 
11996) . For gas densities in galaxy clusters, n, the radiative cool- 
ing time of CR protons is much longer relative to the hadronic 
timescale, Tpp ^ 30 Gyrx(n/10"-' cm"-')"', on which CR protons 
collide inelastically with ambient gas protons of the ICM. This 
is the process in which we are primarily interested here as it pro- 
duces pions that successively decay into synchrotron emitting 
electrons/positrons (so-called secondaries) and gamma-rays. 

There are two principal models that explain radio 
(mini-)halos. In the "hadronic model" the radio emitting 
electrons are produced in ine l astic CR proton - proton (p-p) 
interactions ('Dennison' '1980'; 'Ves trandl 1 198 2*: 'EnBlin et al.l 
1997; Blasi & Colafrancesco 1999; Jpola£& EnBHn 2000; 
Miniatietal. 2001; Miniati 2003; Pfrommer & EnBHn 
2003, 20 04a,,b; Blasi et al. 2007; Pfrommej" et aL 200& 
Pfromm ej 120081: iKushnir et al. 2009; Donne rt et alj l2010alb : 



Keshet & Loeb' 
Fuiita & Ohira 



2010; .Kesheti 
2012) requiring 



only 



,Enfilin et al. 2011 



a very modest frac- 
tion of a few percent of CR-to-thermal pressure. Some clusters, 
known as cool-core clusters, e.g., the Perseus cluster studied 
here, are characterized by a central temperature decrease and a 
correspondingly enhanced central ICM density ('Fabiani ri994l: 
[Hudson et al. 2010). These cool cores provide high target 
densities for hadronic CR interactions, increasing the resulting 
gamma-ray flux. This is the main reason for a larger expected 
flux from Perseus in comparison to oflier nearby clusters 
(iPinzke & Pfrommed l20Tot iPinzke et"an l201 ll) . Another alter- 
native is the "re-acceleration model": during states of powerful 
ICM turbulence, e.g., after a cluster merger, re-acceleration 
might be efficient enough to accelerate CR electrons to high 
enough energi es (~ 10 GeV) that they produce observable 
radio emission (Schlickeiser et al. 1987; Giovanniniet al. 1993", 
^itti et al. 2002; Brunetti & Blasi 2005; Brunetti & Lazarian 
20071 i20TTt lBrunettietal.1 12009|) . This, however, requires a 



around 100 MeV which might be maintained by re-acceleration 
processes by particle-plasma wave interactions at a rate faster 
t han th e cooling processes. We refer the reader to EnBlin et al] 
(1201 ih for a discussion on the strengths and weaknesses of each 
of the models. 

That magnetic fields with central field strengths ranging be- 
tween B ~ {I - 20) fiG are ubiquitous in galaxy clusters and 
permeate the ICM is sugge sted by Faraday rotation measure- 
ments (RM; Kim et al. 1991 ;lClarke et al.l2001l;ICarini & Tavlod 
120021: rv^gt & Enfihni,2005i:lKuchar & EnBlinll201 ih . Combining 
the diffuse synchrotron radio emission with detections of galaxy 
clusters in hard X-rays can also be used to constrain the intra- 
cluster magnetic field if the hard X-ray emission is due to 
inverse Compton up-sc attering of cosmic m icrowave photons 
by CR electrons (see iRephaeh et alj l2008l for a recent re- 
view). The derived magnetic field strengths, typically 0.1 //G, 
are in direct conflict with the much stronger field strengths 
estimated from Faraday RM analyses. However, more sensi- 
tive observations by Swift/BAT have not found evidence of 
a hard power-law tail above the thermal emission. With the 
potential exception of the Bullet Cluster (i.e., IE 0657-558), 
the emission is compatib l e with a multi- temperature plasma 
dAieUo et all l2009l l2010t IWik et all 12011 ). disputing earlier 
claims of the existence of a nonthermal component. Combining 
hadronically-induced synchrotron emission by secondaries from 
CR interactions with the associated pion-decay gamma-ray 
emission would yield significant information on cluster mag- 
netic fields (e.g.. lPfroriirner&EnBUr3l2004btlPfrommeil 120081: 
Ijeltema & Profumdl201lh . Limits on the gamma-ray emission 
provide limits on the intra-cluster magnetic field if the hadronic 
model is applicable to the radio (mini-)halo emission. 

Very high-energy (VHE, £ > 100 GeV) gamma-ray observa- 
tions are very important for delimiting the CR pressure. A sub- 
stantial CR pressure contribution can bias hydrostatic equilib- 
rium cluster masses as well as the total Comptonization parame- 
ter due to the Sunyaev-Zel'dovich effect which i s proportional to 
the thermal energy of a galaxy cluster (e.g., see lPfrommer et aP 
2007). This would at least complicate, if not render it impossi- 
ble, to use these methods for estimating cosmological parame- 
ters with galaxy clusters. There are two ways to assess the bias 
due to nonthermal pressure contributions from CRs, magnetic 
fields or turbulence. At the cluster center a comparison of X- 
ray and optically-inferred gravitational potential profiles yields 
an upper limit of 0.2-0.3 for the ratio of nonthermal pressure to 
the thermal gas pressure dChurazov et alj|2"008ll2010) . On large 
scales, there are indications that hydrostatic mass estimates can 
be biased by up to ~ 20% relative to weak gravitational lens- 
ing masses at R^qo — the radius of a sphere enclosing a mean 
density that is 500 t imes the critical density of the Universe 
dMahdavi et al.ll2008l) . While this could be due to nonthermal 
pressure, weak lensing mass estimates must be interpreted with 
caution since they are intrinsically biased: cluster mass (shear) 
profiles at large radii attain sys tematic modifications due to sub- 
structure and halo triaxiality dBecker & Kravtsovll201 lb . VHE 
gamma-ray emission is a complementary method for constrain- 
ing the pressure contribution of CRs that is most sensitive to 
the cluster core region. Multi-frequency constraints limited the 
CR pressure component to be belo w a few percent of the ther- 



sufficiently long-lived CR ICM electron population at energies 



mal pressure in the best cases (e.g.,|Pf rommer & EnBlinll2004a|; 
iReimer et alj2004l;lAckermann et alTl olOa; Pinzke et al. 201 if. 
However, this approach assumes that the CR component is fully 
mixed with the ICM and is almost insensitive to a two-phase 
structure of CRs and the thermal ICM. In addition to this cos- 
mological motivation, constraining the CR pressure contribu- 
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tion in clusters in combination with cosmological hydrodynam- 
ical simulations provides limits on the average acceleration ef- 
ficiency at strong formation shocks, the CR transport properties 
dEnBlin et alJ 1201 1 ) and the impact of the cluster's dynamical 
state on the CR population. 

Following up on earlier high-energy gamma-ray observa- 
tions of galaxy clusters, in this work we present the deepest 
yet VHE gamma-ray observation of a galaxy clust er (for space- 
based cluster observations in th e GeV-band, see iReimer et alJ 
120031: lAckermann et alJ l2010allbt for ground-based observa- 
tions in the energy band > 10 GeV, see Perkins et al. 200a 
Perkins 2008; Aharonian et al. 2009a b; Domainko et al. 2009^ 
,Galante et al. 2009; Kiuchi et al, ,2009; Acciari et al, 2009^ 
Aleksic et al. 201Qab. We are able to constrain simulation-based 
models for the CR induced gamma-ray emission and the cluster 
magnetic field in the hadronic model for radio (mini-)halos. 

The MAGIC Telescopes observed the Perseus cluster in 
stereoscopic mode from October 2009 to February 2011 for 
a total of ~ 85 hr of effective observation time. This cam- 
paign produced the VHE detection of the head-tail radio galaxy 
IC 310, which lies in the field of view of the Perseus observation 
2010bh and of the central radio galaxy NGC 1275 



(Aleksic et al 



dAleksic et al 



1201 2ah . Here, we focus on the CR induced 



gamma-ray emission from the cluster itself. Galaxy clusters 
are also very promising targets for constraining the da rk mat- 
ter annihilation cross section or decay rate C Colafrances co et alj 
2006HPinzke et alJl2009l:IJeltema et al.ll2009 ■Jcuesta et a l.ll20lU 



Duggeretal.1 l2010t ISanchez-Conde et al.l 1201k iPinzke et al.l 
201 U iGaoet al.ll2012h . We defer any such analysis to future 



work. The presence of a central gamma-ray emitting source 
at energies < 600 GeV, combined with the expected flat dark 
matter annihilation emission profile out to the virial radius 
owing to the substructures that dominate the cluster emis- 
sion (Sanchez-Conde et al. 2011; Pinzke et al. 2011; Ga o et al] 
I2OI2I) . calls for novel analysis techniques that are beyond the 
scope of the present study. 

This work has two companion papers. One is dedicated to 
the NGC 1275 detection using the August 2010 - February 201 1 
campaign data ( Aleks ic et al. 2012a), while the other will ad- 
dress the NGC 1275 multi-wavelength emission from radio to 
VHE and give a comprehensive interpretations of all collected 
data (MAGIC Coll., in prep.). The paper is organized as follows. 
In Sect. 121 we explain and present our target selection and put 
the observation of the cool-core cluster Perseus into the broader 
context of nonthermal cluster radio emission. Specifically, we 
compare Perseus to the expected TeV characteristics of a merg- 
ing cluster, namely Coma. In Sect. [3] we describe the MAGIC 
telescopes, present the observations, and detail the data analy- 
sis and results. These data are used to constrain the CR pop- 
ulation in the Perseus cluster and the magnetic field distribu- 
tion in the hadronic model of the radio mini-halo (Sect. |4|. 
Finally, in Sect. |5] we summarize our findings. All quantities 
are scaled to the currently favored value of Hubble's constant 
Hq = 70kms-i Mpc-i. 



2. Target selection: comparing the nonthermal 
emission of cool-core versus merging clusters 

The Perseus cluster was selected for the MAGIC observa- 
tions as it is the most promising target for the detection of 
gamma-rays coming from neutral pion decay resulting from 
hadronic CR interactions with the ICM dAleksic et al.l l2010at 
iPinzke & Pfrommed I20T0I iPinzke et all 1201 ih . Thisduster of 



galaxies (also known as Abell 426), at a distance of 77.7 Mpc 
(z = 0.018), is the brightest X-ray cluster with a luminos- 
ity in the soft X-ray band ( ranging from 0.1 - 2.4 ke V), 
Lx,o.i-2.4 = 8.3 X lO^^ergs ' dReiprich & Bohringeill2002h . It 
contains a massiv e cool core with high central gas densities of 
about 0.05 cm^^ dChurazov et al.ll2003h . The cluster also con- 
tains a luminous radio mini-halo with an extension of 200 kpc 
(Pedlar et al. 1990). More details on possible VHE gamma-ray 
emission in galaxy clusters can be found in the first MAGIC pa- 
per about the Perseus cluster observation with a single telescope 
(lAleksic et aT]|2010al) . 

One of the most important questions in studies of non- 
thermal cluster emission is the origin of giant radio halos in 
merging clusters and radio mini-halos in cool-core clusters and 
whether they have a common physical origin. Hence it is in- 
structive to compare the prospects for detecting VHE gamma- 
ray emission from two representative clusters of each model 
class, namely Coma and Perseus. First, assuming universality of 
the gamma-ray spectrum. Coma is expected to be fainter than 
Perseus by approximately a factor of 3.4. This i s mainly be- 
cause of the lower central gas density in Coma dPinzke et al.l 
Second, the Coma X-ray emission is more extended than 
in Perseus by the ratio of the half-flux radii of 0.1 8 70.11° - 1.6 
(Pinzke & Pfrommer 2010). Since the sensitivity of Imaging 
Atmospheric Cherenkov Telescopes (lACTs) drops almost lin- 
early with source extension for the angular scales considered 
here, this makes it even more challenging to detect the Coma 
cluster. Third, the CR spectral index is expected to steepen due 
to CR transport processes s uch as stre aming and momentum- 
dependent diff'usion (Enfilin et al.ll201 1I) . The authors of that pa- 
per indicate that these mechanisms alleviate some tension be- 
tween the hadronic model of radio halos and their observations, 
such as Coma. The following line of arguments shows that this 
third point slightly favors Coma as a target source, but the sum 
of all points clearly favors Perseus: 



- Coma should have a central magnetic field strength of 
Bo__^__4//^nc yuG according to Faraday RM studies (e.g., 
iBonafede e t al. 2010). While GHz-radio observations probe 

2.5 GeV electrons that are produced in hadronic interactions 
of 40 GeV CR protons, gamma-rays at 200 GeV are from 

1.6 TeV CRs: hence, for Coma, we have to extrapolate a fac- 
tor of 50 in energy to connect the CRs probed by radio halo 
observations with those probed by lACTs and we expect a 
weak bias due to possible CR spectral steepening. 

- For Perseus, we expect a stronger central magnetic field 
(~ 20/iG, see discussion in Sect. 14.41 1 from adiabatic mag- 
netic field compression during the formation of the cool core. 
While GHz-radio observations probe 1 .2 GeV electrons and 
the 20 GeV CR protons that generate those electrons hadron- 
ically, the MAGIC gamma-ray constraints at 1 TeV corre- 
spond to 8 TeV CR protons (since the detection of VHE 
gamma-ray emission from the central AGN, NGC 1275, ren- 
ders it difficult to observe diffuse emission at lower energies). 
Hence the CRs probed by the radio observations and those 
probed by MAGIC are a factor of 400 different in energy. 
A possible CR spectral steepening due to, e.g., momentum- 
dependent CR diffusion would induce a larger bias in the CR 
pressure and magnetic field constraints. Assuming a change 
in the spectral index by 0.2 between 20 GeV and 8 TeV im- 
plies a decrease of the VHE gamma-ray flux by a factor of 
3.3 — 1.6 times larger than the corresponding decrease of 
flux for Coma (over the eight times smaller energy range). 
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RA [h] RA [h] 

Fig. 1. Left panel: significance skymap of the Perseus cluster above 150 GeV. The positions of NGC 1275 (white cross), which 
coincides with the cluster center, and the head-tail radio galaxy IC 310 (black cross) are marked. Right panel: significance skymap 
of the Perseus cluster region above 630 GeV. The central white circle marks the 0.15° region used to derive constraints (see main 
text for details). The quoted energy thresholds are obtained assuming a spectral index of 4 as observed in NGC 1275. 



In summary, it appears that Perseus is the most promising 
target for detecting CR-induced gamma-ray emission and hence 
we chose this cluster for our deep VHE gamma-ray campaign. 
Assuming that the energy dependence of the CR transport and 
the associated spectral steepening is representative for clusters, 
we would need ten times longer integrations to detect gamma- 
ray emission from Coma in comparison to Perseus. However we 
stress that cool-core and non-cool core clusters are complemen- 
tary and equally deep or deeper observations of merging clusters 
that host giant radio halos — such as Coma — are needed. We 
also emphasize the need to theoretically model how the variation 
of the CR spectral index due to energy dependent CR transport 
processes, e.g., CR diffusion, depends on the mass and dynam- 
ical state of the cluster This will be critical for modeling the 
gamma-ray brightness of clusters and how to interpret the un- 
derlying CR physics. 



3. MAGIC observations, analysis, and results 

The MAGIC telescopes are two 17 m dish lACTs located at the 
Roque de los Muchachos observatory (28.8°N, 17.8°W, 2200 m 
a.s.l.), on the Canary Island of La Palma. Since the end of 2009 
the telescopes are working together in stereoscopic mode which 
ensures a sensitivity of ~ 0.7% of the Crab Nebul a flux above 
approx imately 600 GeV in 50 hr of observations ( Aleksic et al.l 
I2012bl) . For the energies of interest here (i.e., above 600 GeV), 
the point spread function (PSF), defined as a 2-dimensional 
Gaussian, has a cr ^ 0.06°. 

The Perseus cluster region was observed by the MAGIC 
telescopes from October 2009 to February 2011 for a total of 
about 99 hr During the October 2009 - February 2010 campaign 
(45.3 hr), the data were taken in the so called soft-stereo trig- 
ger mode with the first telescope (MAGIC-I) trigger working 
in single mode and the second telescope (MAGIC-II) recording 
only events triggered by both telescopes. The soft-stereo trigger 
mode may result in slig htly deg raded performance at the lowest 
energies with respect to lAleksic et al.l (12012bl) . but has a negli- 
gible impact at the energies of interest here. During the August 
2010 - February 201 1 campaign (53.6 hr), data were taken in the 
standard full-stereo trigger mode (where events are triggered si- 
multaneously by both telescopes). O bservations were performed 
in wobble mode dFomin et al.lll994|) tracking positions 0.4° from 
the cluster center, at zenith angles ranging from 12° to 36°. 



The selected data sample consists of 84.5 hr of effective ob- 
servation time. The data quality check resulted in the rejection 
of about 14.4 hr of data, mainly due to non-optimal atmospheric 
conditions. The standard MAGIC stereo analysis chain was used 
for calibration and image cleaning (see Aleksi c et al.ll2012bl for 
details). In order to combine data taken with different trigger 
modes, we applied a high cut on the shower image size (to- 
tal signal inside the image). Only events with image size above 
150 photo-electrons in both telescopes were kept (the standard 
analysis uses 50 photo-electrons cuts). With this cut, the rate and 
image parameter distributions of the background events are com- 
patible between the two samples. 

The left panel of Fig. [T] shows the significance skymap ob- 
tained with this analysis with an energy threshold of 150 GeV. 
The two AGNs detected during the campaign can clearly be seen 
(NGC 1275 in the center and IC 310 in the lower-right part). 

In this work we limit our attention to pion-decay gamma- 
ray emission resulting from hadronic CR interactions with 
thermal protons of the I CM. Cosmological simulations of 
iPinzke & Pfrommerl (1201 Oh suggest that the spectral energy dis- 
tribution of gamma-rays follows a power-law, F oc £ with 
a spectral index of g - 2.2 si the energies of interest here 
(lAlglyic gt_al.j20 1 Qab ■ The simulated signal is extended with ap- 
proximately 60% of the emission coming from a region cen- 
tered on NGC 1275 with a radius of 0.15° (see Fig. 13 of 
iPinzke & Pfrommeil I20T0I) . The emission from NGC 1275 is 
dominan t below aboiit 6 00 GeV and with a spectral index of 
about 4 (Aleksic et al.l2012a) . Therefore, since the expected CR- 
induced signal is much harder than the measured NGC 1275 
spectrum, we expect it to appear at higher energies with no 
break or cut-off in the energy range covered by MAGIC. Hence, 
we limit the data analysis to e nergies > 630 GeV for which 
the NGC 1275 signal vanishes ( lAleksic et a"l]|2012al) . The right 
panel of Figure[T]shows the significance skymap above 630 GeV. 
In contrast to NGC 1275, the spectrum of IC 310 i s very hard 
and remains detectable above 600 GeV ( lAleksic et al. 2010b|). 
IC 310 is ~ 0.6°(- 10 PSF) away from the cluster center, far 
enough so that its emission does not leak into the signal re- 
gion. However, if not explicitly accounted for, it could affect the 
estimated background. Here, the background is measured with 
three off-source positions at 0.4° from the the camera center and 
> 0.28° away from IC 310. This distance guarantees that there 
is no contamination from IC 310. 
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Table 1. Integral flux ULs. 







N„„ = 325;N^ii = 319.2±10.5 






N„ = 5.8 


t 






Significance (Li&Ma) = 0.28a 




m 








, , 1 .... 1 





0.1 



0.2 



0.3 



0.4 



0.5 



[deg^] 



Fig. 2. &^ distribution above 630 GeV of the signal (data points) 
and background (gray filled region) at the cluster center. The 
dashed line represents the signal region cut within which the ex- 
cess event number is estimated. Here it corresponds to 0.02 deg^. 
The quoted energy threshold is obtained assuming a spectral in- 
dex of 4 as observed in NGC 1275. 



In Fig.|2] we show the 9^ distributioiflof the signal and back- 
ground at the cluster center for energies above 630 GeV. Since 
the emission is expected to be more extended than the MAGIC 
PSF we must optimize the cut in the 6^ distribution which de- 
fines the signal region. We take the following approach. First, 
we estimate the &^ distribution of the expected signal considering 
the simulated surface brightness of the Perseus cluster emission 
smoothed with the MAGIC PSF. Second, we estimate the back- 
ground level from our data. Finally, we calculate the significance 
for signal detection (according to the Eqn. 17 in Li & Mal ll983h 
as a function of the &^ cut. The derived optimal 0- cut was close 
to 0.02 deg^. As shown in Fig.|2] we find no significant excess 
within 0.02 deg^. 

Since we did not detect a signal, we derived integral flux up- 
per limits (ULs) for several energy thresholds. The eff'ective area 
of MAGIC is calculated using point-like source Monte Carlo 
simulation and the ULs have to be corrected to take into account 
the expected source extension. To calculate this correction factor, 
we compare the fraction of the total events inside the signal re- 
gion for a point-like source and for the expected CR induced sig- 
nal. Therefore, the presented ULs can be compared with the the- 
oretical expectations for the region within a radius o f 0.15°. The 
UL estimation is performed using the Rolke method ( iRolke et"an 
|2005) as in Aleksic et al.. ( 201 1) with a confidence level of 95% 
and a total systematic uncertainty of 30%. In Table[l] we present 
the integral flux ULs above specified energy thresholds calcu- 
lated for a spectral index of 2.2. We also give the point-like ULs, 
significance and ULs in number of events. The integral ULs for 
energies above Eth - 1 TeV corresponds to the best sensitivity 
for sources with spectral index 2.2 and it is the most constraining 
value; for this reason we will adopt this UL for the discussion. 
For purposes that will be clear in the following, we also recalcu- 
late the 0.15° integral flux UL above 1 TeV for spectral indexes 
of 2.1,2.3 and 2.5 which are 1.37, 1.38 and 1.39x10"'^ cm^^ s"', 
respectively. 



E,h[GeV] 


" LiMa 




pPL 




630 


0.59° 


84.7 


2.93 


3.22 


1000 


0.15 


41.4 


1.25 


1.38 


1600 


0.33 


38.7 


1.07 


1.18 


2500 


0.38 


28.8 


0.79 


0.87 



Note. Integral flux ULs Ful for a power-law gamma-ray spectrum 
with spectral index 2.2 above a given energy threshold E,i,, both for 
a point-like source (PL) and for a 0.15° extended region, in units 
of 10"'^^ cm"^ s"'. We additionally show the corresponding signifi- 
cance cr^j^^ and ULs in number of events N^^ (before applying the 
source extension correction). " Note that the significance reported for 
E,ii = 630 GeV is slightly different than in figure |2] because different 
cuts were applied. In particular, while computing the 9^ distribution we 
used a hard gamma-ray selection cut, which is normally adopted for de- 
tection purpose, in computing the flux ULs we used softer cuts in order 
to reduce the systematic errors. 



4. Implications for cosmic rays and magnetic fields 

We follow different approaches in constraining the CR pressure 
distribution in the Perseus cluster with the MAGIC ULs. This 
enables us to explore the underlying plasma physics that pro- 
duce the CR distributio n and hence reflec t the Bayesian pri- 
ors in the models (see iPinzke et alj 1201 ll for a discussion). 
Our approaches include (1) a simplified analytical approach 
that assumes a constant CR-to-thermal pressure and a momen- 
tum power-law spectrum (Sect. 14.11 ). (2) an analytic model of 
CRs derived from cosmological hydrodynamical simulations of 
the formation of galaxy clusters (Sect. 14.21 ). and (3) the use of 
the observed luminosity and surface brightness profile of the 
radio mini-halo in Perseus to place a lower limit on the ex- 
pected gamma-ray flux in the hadronic model of radio mini- 
halos — where the radio-emitting electrons are secondaries 
from CR interactions (Sect. 14.3b . This provides a minimum CR 
pressure which, using tight gamma-ray limits/detections, checks 
the hadronic model of the formation of radio mini-halos. (4) 
Alternatively, by constructing a CR distribution that is just al- 
lowed by the flux ULs, and requiring the model to match the ob- 
served radio mini-halo data, we can derive a lower limit on the 
magnetic field strength (Sect. 144b . Note, however, that this limit 
assumes that the observed synchrotron emission is produced by 
secondary electrons resulting from hadronic CR interactions. 



4.1. Simplified analytical CR model 

Here, we adopt a simplified analytical model that assumes a 
power-law CR momentum spectrum, / oc p^", and a constant 
CR-to-thermal pressure ratio, i.e., we adopt the isobaric model 
of CRs following the approach of IPfrommer & EnBlinI (l2004ai) . 
To facilitate comparison with earlier analytical work in the lit- 
erature, we do not impose a low-momentum cutoff on the CR 
distribution function, i.e., we adopt q - Q (this assumption 
can b e easily generalized to an arbitrary q using e.g.. Fig. 1 of 
Pfr ommer & EnBlinI l2004l . Since a priori the CR spectral in- 
a, is unconstrained, we vary it within a plausible range of 
2.1 < a < 2.5. The central value of this spectral range is compat- 
ible with the radio spectral index in the core of the Perseus radio 



' The is the squared angular distance between the arrival directio n 
of events and the nominal source positions (see e.g. lDaum et alJ[T997h . 



- The hadronic interaction physics guarantees that the CR spectral 
index coincides with that from pion-decay gamma-ray emission. 
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Table 2. Constraints on CR-to-thermal pressure ratio. 



a 


^CR.max [%] 


^CR,rain [%] 


-^CR.max/^CR.min 


F" - 

y.min 


2.1 


0.77 


0.42 


1.8 


7.4 


2.2 


1.12 


0.35 


3.2 


4.3 


2.3 


2.17 


0.38 


5.7 


2.4 


2.5 


11.6 


0.67 


17.3 


0.8 



Note. Constraints on CR-to-thermal pressure ratio in the Perseus cluster 
core, XcR.max> which is assumed to be constant in the simplified analyti- 
cal model. Those constraints are compared to minimum CR-to-thermal 
pressure ratios, XcR.min, and minimum gamma-ray fluxes in the hadronic 
model of the Perseus radio mini-halo. ° Minimum gamma-ray flux in the 
hadronic model of the Perseus radio mini-halo, F^^in(> 1 TeV), in units 



mini- haloof Q-v = a/2 ~ 1.25 (ISiibringlll99 3')PI assuming a cen- 
tral magnetic field strength of 20 fiG, appropriate for cool-core 
clusters, the CR protons responsible for the GHz radio emitting 
electrons have an energy of ~ 20 GeV and are ~ 400 times less 
energetic than the CR protons with an energy of 8 TeV that are 
responsible for gamma-ray emission at 1 TeV. This is consistent 
with the concavity in the universal CR spectrum that connects 
the steep-spectrum low-energy CR population around GeV en- 
ergies (a ^ 2.5) with the harder CR population at TeV ener- 
gies (a ^ 2.2) (Pinzke & Pfrommer 2010). To model the ther- 
mal pressure, we adopt the measured el ectron temperature an d 
density profiles for the Perseus cluster jChurazov et alj|2003h . 
The temperature profile has a dip in the central cool core region 
and otherwise a constant temperature of kT - 1 keV. The den- 
sity profile is hybrid in the sense that it combines Einstein X-ray 
observations on large scales wit h high-resolution XM M-Newton 
observations of the cluster core ( Churazov et alj|2003i) . 

Table|2]shows the resulting constraints on the CR-to-thermal 
pressure ratio, Xqr - (Pcr) I (Ah), averaged within the virial 
radius, 7?vii = 2 Mpc, that we define as the radius of a sphere 
enclosing a mean density that is 200 times the critical density of 
the Universe. The inferred constraints on Xqr strongly depend 
on a due to the comparably large lever arm from GeV-CR en- 
ergies (that dominate the CR pressure) to CR energies at 8 TeV. 
Using the integral flux UL above 1 TeV, we constrain Xcs. to 
be between 0.71% and 11.6% (for a varying between 2.1 and 
2.5). For a spectral index of 2.2, favored by the simulation-based 
model of Pinzke & Pfrommer, (i2010i) at energies > 1 TeV, we 
obtain XcR < 1.1%. 

4.2. Simulation-inspired CR model 

For a more realistic approach, we turn to cosmological hydrody- 
namical simulations. These have considerable predictive power, 
e.g., they calculate the CR spectrum self-consistently rather than 
leaving it as a free parameter, and permit tests of various assump- 
tions about the underlying CR physics. We adopt the universal 



' The reported radi al spectral steepening of the radio mini-halo emis- 
sion jGitti et alj2002) could be an observational artifact owing to a poor 
signal-to-noise ratio in the outer core of the cluster and the ambiguity in 
determining the large scale Fourier components owing to nonuniform 
coverage of the Fourier plane and missing short-baseline information: 
the so-called "missing zero spacing"-problem of interferometric radio 
observations. By comparing the spectral index distribut ion of the three 
radio maps (at 92 cm, 49 cm, and 21 cm, ISiibrinsI 1 993h . radial spectral 
flattening depending on the chosen radial direction seems possible. 



spectr al and spatial CR model developed bv lPinzke & Pfrommer! 

that uses only a density profile as input which can be 
inferred from cosmological simulations or X-ray observations. 
This analytic approach models the CR distribution, and the as- 
sociated radiative emission processes from hadronic interac- 
tions with gas protons, from radio to the gamma-ray band. The 
gamma-ray emission from decaying neutral pions dominates 
over the IC emission from prim ary and sec ondary electrons 
above 100 MeV in clusters (e.g., lPfrommeii f2008). This ana- 
lytic formalism was derived from high-resolution simulations 
of clusters that included radiative hydrodynamics, star forma- 
tion, supernova feedback, and followed CR physics spectrally 
and spatially by tracing the most important injection and loss 
processes self-consistently while accounting for the CR pressure 
in the equation of motion dPfrommer et al.l 120061 : lEnBlin et al.l 
2007; Jubelgas et al. 20Q3). The results are consistent with ear- 
lier numerical results within their range of applicability regard- 
ing the overall characteristics of the C R distributi on and the 
associated radiative emission processes (iDolag & En filin 200d; 
Miniati et al !''2001l; iMiniatH 120031: iPfrommer et all 12007. qool 
Pfrommer 2008). 

There are (at least) two major uncertainties in modeling 
the CR physics that significantly affect the resulting spatial and 
spectral CR distribution (and may impact the universality of the 
CR populations across different cluster masses): the CR accel- 
eration efficiency and the microscopic CR transport relative to 
the thermal plasma. We will discuss each process separately. 
The overall normalization of the CR and gamma-ray distribu- 
tion scales with the maximum acceleration efficiency at structure 
formation shocks. Following recent observations (Helder et al.l 
2009) and theoretical studies ( Ka ng & Jones.2005i) of supernova 
remnants, we adopt an optimistic but realistic value of this pa- 
rameter and assume that 50% of the dissipated energy at strong 
shocks is injected into CRs while this efficiency rapidly de- 
creases for weaker shocks. The vast majority of internal for- 
mation shocks (merger and flow shock s ) are weak with Mach 
numbers M < 3 (e.g., iRvu et"al] 120031: pfrommer etal] 120061) 
which we presume have low CR acceleration efficiencies and 
do not contribute significantly to the CR population in clusters. 
Instead, strong shocks during the cluster formation epoch and 
strong accretion shocks at the present time (at the boundary of 
voids and filaments/super-cluster regions) dominate the injection 
of CRs which are adiabatically transported throughout the clus- 
ter. These models provide a plausible UL for the CR contribution 
from structure formation shocks in galaxy clusters which can be 
scaled with the effective acceleration efficiency. Other possible 
CR sources, such as starburst-driven galactic winds and AGNs, 
could increase the expected gamma-ray yield, but have been ne- 
glectedQ 



* Adiabatic CR losses during the expansion phase of the galactic 
winds imply a negligible contribution to the CR energy density of 
the ICM today (Aleksic et al. 2010a) but integrated over cluster his- 
tory, they may contribute an interesting seed population that however 
needs to be re-accelerated. It is currently unclear whether AGN jets 
are p owered by the Poynt ing flux or by the ki netic energy of hadron s 
fe.g.. lCelotti et al. 1998: H irotani et a l. 1998: Si k"ora"& Madeiskil2000h . 
Even if there was a non-negligible CR ion component, it would 
probably be confined to dilute radio lobes that rise almost unper- 
turbed in the cluster atmosphere to the less dense outer cluster re- 
gions tMcNam ara & Nulsen 2007) where these lobes may be destroyed 
by shocks and dissolve d within the ICM ( EnBlin & Briiggen 20Q2; 
IPfrommer & Jonesll201 ih . Since the CRs are then released into an am- 
bient ICM with a small density, this implies a negligible contribution to 
the cluster gamma-ray emission due to pion decay. 
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These cosmological simulations only consider advective 
transport of CRs by turbulent gas motions which produces a cen- 
trally enhanced profile. However, active CR transport such as CR 
diffusion and streaming flattens the CR radial profile, produc- 
ing a spatially constant CR number density in the limiting case. 
Although the CR streaming velocity is larger than typical advec- 
tion velocities in a galaxy cluster, it becomes comparable to or 
lower than this only for periods with trans- and supersonic clus- 
ter turbulence during a cluster merger. As a consequence a bi- 
modality of the CR spatial distribution results with merging (re- 
laxed) clusters s howing a centrall y concentrated (flat) CR energy 
density profile (lEnBlin et al.ll201 1 ). This produces a bimodality 
of the diffuse radio and gamma-ray emission of clusters, since 
more centrally concentrated CRs will find higher target d ensities 
for hadronic CR proton interactions (lEnBlin et al.ll201 11) . Thus, 
relaxed clusters could have a reduced gamma-ray luminosity by 
up to a factor of five. The goal of this and future searches for 
gamma-ray emission from clusters is to constrain the accelera- 
tion physics and transport properties of CRs. 

In the following, we use two different CR models. In our 
optimistic CR model (simulation-based analytics with galax- 
ies), we calculated the cluster total gamma-ray ffux within a 
given solid angle. In contrast, we cut the emission from indi- 
vidual galaxies and compact galactic-sized objects in our more 
conservative model (simulation-based analytics without galax- 
ies). In short, the ICM is a multiphase medium consisting of 
a hot phase which attained its entropy through structure for- 
mation shock waves dissipating gravitational energy associated 
with hierarchical clustering into thermal energy. The dense, cold 
phase consists of the true interstellar medium within galaxies 
(which is modelled with an effective equation of state following 
ISpringel & Hernquistll2003b and at the cluster center as well as 
the ram-pressure stripped interstellar medium. These cold dense 
gas clumps dissociate incompletely in the ICM due to insuffi- 
cient numerical resolution as well as so far incompletely un- 
derstood physical properties of the cluster plasma. All of these 
phases contribute to the gamma-ray emission from a cluster To 
assess the bias associated with this issue, we performed our anal- 
ysis with both limiting cases bracketing the realistic case. 

We again adopt the hybrid electron density and temper- 
ature profile from X-ray observations of the Perseus cluster 
dChurazov et al. 20030, modifying the temperature profile in the 
outer cluster regions beyond 0.2/? 200 - 400 kpc. This produces 
the characteristic decline toward the cluster periphery in accord 
with a fit to the u niversal profile derived from cosmological 
cluste r simulations (iPinzke & Pfronuneill2010t iPfrommer et al.l 
l2007h that also follows the behavior of a nearby sample of deep 



^ While the X-ray data show a complicated ICM structure including 
cavities blown by the central AGN, the use of deprojected, spherically 
averaged gas density and temperature profiles can be justified as fol- 
lows. The inhomogeneities of the ICM can be quantified statistically 
using a clumping factor Cp = -^(p^)/ (p) that estimates the relative 
bias factor to true density profile. However, since we are interested 
in the CR-induced gamma-ray emissivity which also scales with the 
square of the gas density (modulo a dimensionless CR-to-gas ratio that 
we quantify by using CR simulations) our modeling should inherit ex- 
actly the same bias as in X-ray-inferred density profiles. We emphasize 
that the overall X-ray emission is smooth with fluctuations on the level 
< 20 - 25% (Fabia n et alj|2011 ). Those substructures are only ampli- 
fied by unsharp mask imaging. Moreover, the X-ray cavities (interpreted 
as AGN bubbles) do not correlate with structures in the radio mini-halo 
and only make up a small fraction of the radio-emitting volume (< 10"^ , 
since the y only subtend a li near scale of < 10% relative to the radio 
emission. iFabian et alj20H"h . justifying our procedure. 



Perseus sim. -based analytics w/ gal. 

sim. -based analytics w/o gal. 

T sim. -based analytics w/o gal. (xO.8) 




1000 
E [GeV] 



Fig. 3. Integral ffux ULs of t he single telescope o bservations 
(point-like ULs, upper arrows; I Aleksic et al.ll2010ah and of this 
work (solid arrows; Table 1) are shown. We compare ULs with 
the simulated integrated spectra of the gamma-ray emission from 
decaying neutral pions that result from hadronic CR interactions 
with the ICM in the Perseus cluster coming from within a ra- 
dius of 0.15° around the center The conservative model with- 
out galaxies (solid line) is shown together with the model with 
galaxies (dotted line). The ffux UL at 1 TeV is a factor of 1.25 
below the conservative model and hence constrains a combina- 
tion of CR shock acceleration efficiency and CR transport pro- 
cesses. We also show the minimum gamma-ray ffux estimates 
for the hadronic model of the Perseus radio mini-halo (dashed 
line with minimum ffux aiTows) using the universal gamma-ray 
spectrum resulting from pion decay (Pinzke & Pfrommer 2010|) 
and adopting the spectral index, 2.2, of this spectrum around TeV 
energies. 

Chandra cluster data (" Vikhlinin et al.ll2005h . WTiile this modifi- 
cation has little inffuence on the expected gamma-ray emission 
(in projection onto the core region) as the densities drop con- 
siderably in these regions, the resulting profiles for the CR-to- 
thermal pressure ratio are changed as we will discuss below. We 
estimate t he spatial and spectr al distribution of CRs using the 
model by iPinzke & Pfromme^ (I2OIO ). Figure [3] shows the ex- 
pected spectrum for Perseus within an aperture of radius 0.15°. 
The MAGIC limit for Ey > 1 TeV falls below the ffux level of 
the conservative model (that does not take into account the emis- 
sion from galaxie^ by 20%; thereby constraining assumptions 
about the adopted CR physics in the simulations and the result- 
ing CR pressure. 

Figure |4] shows the CR-to-thermal pressure ratio, Xqr = 
(^'cr) / (Ah) as a function of radial distance from the Perseus 
cluster center, in units of the virial radius, /?vii. To compute 



* Note that the models "with galaxies" and "without galaxies" dif- 
fer from those used in t he first MAGIC paper on Perseus observations 
with a single telescope i Aleksic et alJl20103) . Previously, the models 
were based on the simulated cool core cluster g51 (which is similar 
in morphology to Perseus) while here we use the universal CR model 
— which is based on high-resolution simulations of 14 galaxy clusters 
that span almost two decades in mass — to compute the expected pion 
decay emission (Pinzke & Pfrommer 2010). The spectral normalization 
(between the two models) therefore increased from 1.5 to 2.4 and the 
gamma-ray flux in the model "without galaxies" used here differs by 
10%. 
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Zcr(«), scaled to 0.15° UL" 
X^^(<R), scaled to 0. 1 5° UL - 
^cr(^) rnodel 
Xf~J<R) model 




Fig. 4. CR-to-thermal pressure ratio, Xqr - (Pcr) I (Ah) at the 
radial distance from the Perseus cluster center, R/Rvh (red lines), 
and integrated up to R/Ryn (bl ue lines), using the simul ation- 
based analytical model of CRs dPinzke & Pfromnie^l2010l) . The 
simulation model (dashed) is contrasted to a model that has been 
scaled by the MAGIC constraints obtained in this work (solid). 



Pcr, we assume a low-momentum cutoff of the CR distribu- 
tion of q = 0.8 nipC found in cosmological simulations by 
iPinzke & Pfromm er (2010). This cutoff derives from the high 
cooling rates for low energy CRs due to Coulomb interactions 
with the thermal plasma. The ratio Xcr rises toward the outer 
regions because of the higher CR acceleration efficiency in the 
peripheral strong accretion shocks compared to the weak cen- 
tral flow shocks. Adiabatic compression of a mixture of CRs 
and thermal gas reduces the CR pressure relative to the ther- 
mal pressure due to the softer CR equation of state. The strong 
increase of Xcr toward the core is a remnant of the formation 
of the cool core in Perseus. During this transition, the mixed 
gas of CRs and thermal particles has been adiabatically com- 
pressed. While the thermal gas radiates on a comparatively fast 
thermal bremsstrahlung timescale, the long hadronic interaction 
time scale for energetic CRs ensures an accumulation of this 
population, thus diminishing the thermal pressure support rel- 
ative to that in CRs. 

The MAGIC flux limits constrain Xqr of the simulation- 
based analytical model to be less than 1.6% within 0.15° ( 200 
kpc). Assuming this spatial CR profile yields a CR-to-thermal 
pressure ratio < 1.7% within /?vii - 2 Mpc and < 5% within 
20 kpc. When scaling the Xcr profile to the flux ULs, we used 
the conservative model of iPinzke & Pfrommeii ( 1201 Oh that ne- 
glects the contribution of cluster galaxies since those increase 
only the gamma-ray yield but do not contribute to the CR pres- 
sure within the ICM. The XcR-limit within the virial radius is 
larger by a factor of 1.5 than a pure power-law spectrum (see 
Sect. 14. II ) because the concave curvature of the simulated spec- 
trum accumulates additional pressure toward GeV energies rela- 
tive to a pure power-law. 

With these gamma-ray flux ULs, we constrain the CR 
physics in galaxy clusters: this either limits the maximum ac- 
celeration efficiency of CRs at strong structure formation shocks 
to < 50% or indicates possible CR streaming and diffusion out 
of the cluster core region. The latter populates the peripheral 
cluster regions and Xcr increases toward the clust er periphery 
at the expense of a decrease of the central Xcr (lEnBlin et alJ 



1201 ih compared to the simulation model. The X-ray morphol- 
ogy of the central region in Perseus shows spiral structure in the 
density and tem perature m aps with an anti-correlation of both 
quantities dFabian et alJl201 1 ). This resembles sloshing motions 
after a past merger event, suggesting that Perseus is currently 
relaxing. This interpretation is supported by rec ent magneto- 
hydrodynamical simulations (IZuHone et al.ll20l T) which argue 
that magnetic fields are dra ped at th e contact discontinuity of 
the remnant cool co re (Lvutikov 2006; Dursi & Pfrommer 200^ 
iPfrommer & Durs^i201 0) that is sloshing in the potential of the 
parent halo. The abrupt drop in the radi o surface brightnes s map 
outside the spiral pattern (see Fig. 8 in lFabian et alj|20Tl ) may 
then be produced by a strong decrease of magnetic field strengths 
outside the sloshing core. If CR streaming and diffusion out 
of the central core region is indeed correlated with a dynami- 
cal relaxation of a cluster after a merger event as suggested by 
lEnBlin et al.l (1201 lb . this would render CR transport more plau- 
sible for explaining the smaller gamma-ray flux relative to the 
simulation model. By the same token, it may argue for a more 
extended gamma-ray emission signature than that seen in the ra- 
dio, further justifying the larger source extension, 0.15°, that we 
adopted in this work (which is twice that of the largest radial 
extend of the mini-halo of 0.075°). 



4.3. Minimum gamma-ray flux 

For clusters with radio (mini-)halos we can derive a minimum 
gamma-ray flux in the hadronic model of radio (mini-)halos. A 
stationary distribution of CR electrons loses all its energy to syn- 
chrotron radiation for strong magnetic fields (B » Bcmb - 
3.2juG, where Bcmb is the equivalent magnetic field strength 
of the cosmic microwave background (CMB) so that R^y^^l%n 
equals the CMB energy density). Hence the ratio of gamma-ray 
to synchrotron flux becomes independent of the spa tia l distribu- 
tion of CRs and thermal gas (Volk 1989; Pohl 1994: IPfrommeil 
2008.) . particularly with the observed synchrotron spectral index 
ay = a 12 ^ 1. This can be easily seen by considering the pion- 
decay induced gamma-ray luminosity, Ly, and the synchrotron 
luminosity, Ly, of a steady state distribution of CR electrons that 
has been generated by hadronic CR interactions, 



Ly - Ay J dV ncRng^^, 



I 



liV«CR«gas 



^(a-v + l)/2 
ecMB + SB 



^ Ay JdVncRngi^s foreB»ecMB- 



(1) 

(2) 
(3) 



Here, sb and scmb denote the energy density of the magnetic 
field and the CMB, respectively. Ay and Ay are dimensional con- 
stants that depend on the hadronic p hysics of the interaction 
(iPfrommer et al.l 120081 ; iPfrommeii 12008) . Hence we can derive 
a minimum gamma-ray flux in the hadronic model: 



y.min 



Ay AnD\ 



(4) 



lum 



where Ly is the observed luminosity of the radio mini-halo and 
D\am - 77 Mpc denotes the luminosity distance to Perseus. 
Lowering the magnetic field strength would require increasing 
the CR electron energy density to reproduce the observed syn- 
chrotron luminosity and thus increases the associated gamma- 
ray flux. The maximum emission radius of the radio mini-halo 
of 100 kpc corresponds to an angular size of 0.075°, well within 
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the MAGIC PSF; hence Ly does not need to be cut to match the 
angular region tested by the MAGIC ULs. 

Hence in this section, we assume B » Bcmb everywhere 
in the radio-emitting region, and (in order to obtain this min- 
imum gamma-ray flux) require a considerable drop in the CR 
distribution outside the radio (mini-)halo. While we do not make 
any assumptions on how to realize such a scenario, we discuss 
possible physical processes that could be causing it. The large 
magnetic field can be realized by adiabatically compressing the 
magnetic field during the formation of the cool core. The edge of 
the radio emitting region can be caused by a region of predom- 
inantly toroidal magnetic field which efficiently confines CRs 
to the central region. Such a magnetic field configurat ion could 
be arr anged by magnetic draping at the sloshing co re ( Lyutikovl 
l2006tlDursi & Pfrommeil2008t (Pfrommer & Dursill201d) or the 
saturated (non-Unear) state of the heat- flux driven buoyancy in- 
stability (Ouataert 2008; Parrish & Ouataert 2008"). The drop in 
CR number density outside this region can be realized by CR 
streaming to the peripheral cluster regions with substantially 
lower target densities that would imply a negligible contribution 
to the gamma-ray flux outside the radio mini-halo (EnBlin et al] 
1201 li) . The next generation of cluster simulations that couple 
CR physics self-consistently to magnetic fields and account for 
anisotropic transport processes of the thermal gas as well as CRs 
will be needed to scrutinize such a presented scenario. 

In contrast to our previous analysis in Aleksic et"al] (1201 Oab . 
we derive here the absolute minimum gamma-ray flux rather 
than the "physical" minimum gamma-ray flux that depends on 
the assumed magnetic field distribution and requires an uncer- 
tain extrapolation into the outer regions of the cluster The den- 
sity and pressure profile of Perseus as derived from X-ray ob- 
servations can both be described by double-beta profiles with an 
inner cor e radius of 57 and 47 kpc, respectively (IChurazov et al.l 
[2003; Pf rommer et al.l l2005). Hence, there is only a small inter- 
val between ~ 50 and ~ 100 kpc, the maximum radius of the 
radio mini-halo, that is sensitive to the outer CR and magnetic 
field profile. Only the extrapolation of that behavior well beyond 
the radio mini-halo to the virial radius will determine whether 
the energy density in the magnetic field becomes implausibly 
large in comparison to that of the thermal gas, thereby possibly 
violating the energy conditions. 

The results for the minimum gamma-ray flux Fymin(> 
1 TeV) and the minimum CR-to-thermal pressure ratio, 
^CR,min = ^CR^'y.min/^y.i.so are presented in Tableland Fig. 
[3] (assuming a - 2.2). Here, Fy iso is the gamma-ray flux in the 
isobaric model of CRs introduced in Sect. 14.11 The ratio of the 
maximum to minimum CR pressures, XcR.max/^^CR.min, varies be- 
tween 1.8 and 17.3 for a spectral index between 2.1 < a < 2.5. 
For the spectral index a - 2.2 of the universal gamma-ray spec- 
trum around TeV energies, the ratio is X^cR.maxZ-'^CR.min = 3.2. 
This puts the long-sought gamma-ray detection of clusters, in 
particular for Perseus, within the reach of deeper campaigns with 
the possibility of scrutinizing the hadronic emission model of ra- 
dio (mini-)halos. 

4.4. Implications for the cluster magnetic field 

As shown in Sect. 14.31 an absolute lower limit on the hadronic 
model gamma-ray emission comes from assuming high mag- 
netic field strengths, B » Bqmb, everywhere in the radio- 
emitting region. Using our UL on the gamma-ray emission (and 
on Pgr) we can turn the argument in Sect. 14.31 around to de- 
rive a lower limit on the magnetic field strength needed to ex- 
plain the observed diffuse radio emission within the hadronic 



model (iPfrommer & EnBUnll2004bl) . Lowering the gamma-ray 
limit will tighten (increase) the magnetic field limit. If this con- 
flicts with magnetic field measurements by means of other meth- 
ods, e.g., Faraday RM, this would challenge the hadronic model 
of radio (mini-)halos. The method that we use to constrain the 
magnetic field, B, depends on the assumed spatial structure of B 
that we parametrize by 



B(r). 



(5) 



as suggested by a cosmological MHD simulation of a cool core 
cluster (Dubois & Tevssier 2008), which demonstrates a tight 
coiTelation of the magnetic field with the ICM gas density and 
highlights the importance of cooling processes in amplifying the 
magnetic field in the core of galaxy clusters up to one order of 
magnitude above the typical amplification obtained for a purely 
adiabatic evolution. Mor eover, such a p arametrization is favored 
by Faraday RM studies dBonafede et al.ll2010l;lKuchar & EnBlinI 
i2011 , and references therein). 

Generally, Faraday RM determinations of the magnetic field 
strength, for instance using background sources seen through 
clusters, depend intrinsically on the magnetic correlation length. 
Recent Faraday RM studies yield estimates for t he central mag- 
netic field of typically 5 //G for merging clusters dBonafede et al.l 
2010, for Coma) and si gnificantly higher value s in cool core 
clusters of around 16 yuG ('Kucha r & EnBlin F201 1, for Hydra A). 
For the Perse us cluster, Faraday RMs are available only on very 
small scales dTavlor et alj|2006l) . i.e., a few tens of pc. The RM 
estimates are ~ 7000 rad m"^ implying magnetic field strengths 
of ~ 25 ^G. This assumes that the Faraday screen is situated 
within the ICM. This location is not, however, certain sin ce vari- 
ations of 10% in the RM are observed on pc-scales C Taylor et al.l 
l2002h . while ICM magnetic fields are expected to be ordere d on 
scales of a few kpc dTavlor et alJl2006l;IVogt & EnBUnll2005h . 

To proceed, we derived a deprojected radio emissivity pro- 
file. We fit the point-source subtracte d, azimuthally ave raged 
surface brightness profile at 1.38 GHz dPedlar et al.ll 19901) with 
a y6-model. 



l+l^ 



-3/3+1/2 



(6) 



where = 2.3 x 10"' Jy arcmin"^, - 30 kpc, and p - 
0.55. This profile is valid within a maximum emission radius 
of 100 kpc. An Abel integral deprojection then provides the ra- 
dio e missivity distribution (see Appendix of Pfrommer & EnBlinI 
I2004ah . 

where S denotes the beta function. To constrain the magnetic 
field, we proceed as follows: 

1 . Given a model for the magnetic field characterized by as and 
an initial guess for Bq, we determine the profile of XcR(r) 
such that the hadronically produced synchrotron emission 
matches the observed radio mini-halo emission over the en- 
tire extent. 

2. For this XcRir) profile, we compute the pion-decay gamma- 
ray surface brightness profile, integrate the flux within a ra- 
dius of 0.15°, and scale the CR profile o match the corre- 
sponding MAGIC UL. This scaling factor, XcRfi, depends on 
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Table 3. Constraints on magnetic fields in the hadronic model. 





Minimum magnetic field, fio.min \pG]: 








a 




2.1 


2.2 


2.3 2.5 


0.3 


5.86 


4.09 


3.15 2.06 


0.5 


8.62 


6.02 


4.63 3.05 


0.7 


13.1 


9.16 


7.08 4.68 




Corresponding 


;XcR(100kpc) [%]: 


0.3 


1.7 


2.5 


4.9 26.7 


0.5 


1.7 


2.5 


4.8 26.1 


0.7 


1.6 


2.3 


4.5 23.6 



Note. Constraints on magnetic fields in the hadronic model of the 
Perseus radio mini-halo and the corresponding CR-to-thermal pressure 
ratio (at the largest emission radius of 100 kpc) such that the model 
reproduces the observed radio surface brightness profile. 

the CR spectral index, a, (assuming a power-law CR popula- 
tion for simplicity), the radial decline of the magnetic field, 
as, and the initial guess for Bq. 

3. We then solve for Bq requiring that it matches the observed 
synchrotron profile while fixing the profile of XcR(r) deter- 
mined in the previous two steps. Note that for Bq » Bcmb 
and a radio spectral index of ffy = 1, the solution is degener- 
ate, as can be seen from Eqn.|2l 

4. Inverse Compton cooling of CR electrons on CMB photons 
introduces a scale that depends on the CMB -equivalent mag- 
netic field strength, Bcmb - 3.2 yuG; if magnetic field values 
are close to this value, the resulting synchrotron emission de- 
pends non-linearly on the inferred magnetic field strengths. 
Hence, we iterate, repeating the previous steps until conver- 
gence for the minimum magnetic field Bq. 

Table [3] gives the resulting lower limits for Bo that depend 
sensitively on the assumptions of a and ob- Its behavior can be 
understood as follows. (1) The hardest CR spectral indices cor- 
respond to the tightest limits on Bq. This is because for an UL 
for CR energies of around 8 TeV (as probed by 1 TeV gamma- 
rays) and a CR population with a softer spectral index (larger a) 
there is a comparably larger fraction of CRs at 25 GeV available 
which produce more radio-emitting electrons. Thus, lower mag- 
netic field strengths can be accommodated while still matching 
the observed synchrotron flux. (2) For a steeper magnetic de- 
cline (larger ag), the CR number density must be higher to match 
the observed radio emission profiles. This would yield a higher 
gamma-ray flux so the ULs are more constraining. This implies 
tighter lower limits for Bq. 

The inferred values for the minimum magnetic field 
strengths in Table[3]range from 2-13 fiG for the values of a and 
Ob used in this study and suggested by radio observations. These 
are much lower than the thermal equipartition value in the cen- 
ter of Perseus, Bgq.o - 80 fiG or magnet ic field estimates f rom 
the Faraday RM values on small scales dTavlor et al.ll2006l) . In 
Table [3] we additionally give the corresponding values for the 
CR-to-thermal pressure ratio (at the largest emission radius at 
100 kpc) such that the model reproduces the observed radio sur- 
face brightness profiled Since they are derived from flux upper 

' Note that in this section, radial profiles of Xcr are uniquely deter- 
mined by the adopted model for the magnetic field and the observed 
synchrotron surface brightness profile of the radio mini-halo. This dif- 



limits they are also ULs on the CR-to-thermal pressure ratio. 
The corresponding values for Xqr in the cluster center are lower 
than 5% for the entire parameter space probed in this study. We 
conclude that there is still considerable leeway for the hadronic 
model as an explanation of the radio mini-halo emission. 

5. Conclusions 

MAGIC observed the Perseus cluster for a total of ~ 85 hr of 
high quality data between October 2009 and February 2011 . This 
campaign res ulted in the detection of the head-tail radio galaxy 
IC 310 (Aleks ic et 10 120 10b) an d the detect ion of the central ra- 
dio galaxy NGC 1275 iAleksic et al.ll20"T2a ). No significant ex- 
cess of gamma-rays was detected from the cluster central region 
at energies above 630 GeV where the NGC 1275 emission van- 
ishes. The flux UL for the CR-induced emission above 1 TeV, for 
a region of radius of 0. 15° around the cluster center, corresponds 
to 1.38 X 10"'^ cm-2 s-'. 

We constrain the CR population and the magnetic field 
strength in Perseus using contrasting models that differ in their 
assumptions and are constructed to circumvent our lack of 
knowledge of the true distributions. (1) We use a simplified "iso- 
baric CR model" with constant CR-to-thermal pressure fraction 
and power-law momentum spectrum — a model that has been 
widely used in the literature. (2) This is complemented by an an- 
alytical approach based on cosmological hydrodynamical sim- 
ulations (Pinzke & Pfrommer 2010) that provide a CR distribu- 
tion, combined with the observed density profile. (3) We finally 
use a pragmatic approach that models the CR and magnetic field 
distributions to reproduce the observed emission profile of the 
Perseus radio mini-halo. 

Within the simplified analytical approach we can constrain 
the CR-to-thermal pressure, Xqr < 0.8% and 12% (for a CR 
or gamma-ray spectral index, a, varying between 2.1 and 2.5). 
For the spectral index at TeV energies of a = 2.2, favored by 
simulations, we find that Xcr < 1.1%. The simulation-based ap- 
proach gives Xcr < 1.7%. This latter value is a factor of 1.5 
less constraining because of the concave curvature of the simu- 
lated spectrum that has higher partial pressure toward GeV en- 
ergies relative to a pure power-law spectrum. This constraint is 
a factor of 1 .25 below the simulation model and - for the first 
time - limiting the underlying physics of the simulation. This 
could either indicate that the maximum acceleration efficiency 
of CRs relative to the total dissipated energy at strong structure 
formation shocks is < 50% (i.e., smaller than the value assumed 
in the simulations) or may point to CR streaming and diffusion 
out of the cluster core region that lowers the central Xcr val- 
ues (lE nBhn et al.ll201 iT) . The observed spiral X-ray structure in 
the central cluster region suggests that Perseus is currently in 
a relaxation state following a past merging event. If a net out- 
ward CR transport is indeed correlated with a dynamical relax- 
ation state of a cluster, this would render CR transport a plausi- 
ble agent that lowers the gamma-ray flux in comparison to the 
simulation model that neglects this mechanism. 

Our third, pragmatic, approach employed two different as- 
sumptions. First, adopting a strong magnetic field everywhere 



fers from the simplified analytical CR model where Xcr = const. 
(Sect l4. l| l and contrasts with the simulation-based model where XcR(r) 
is derived from cosmological cluster simulations (Sect l4.2l l. For most 
of the cases studied here, the CR-to-thermal pressure ratio profiles vary 
within a factor of two, and flat profiles of both the relative magnetic and 
CR pressure profiles are able t o reproduce the observed ra dio surface 
brightness profile (see Fig. 2 of lPfrommer & Enfilinll2004bh . 



10 



Aleksic et al.: Constraining Cosmic Rays in the Perseus Cluster with MAGIC 



in the radio-emitting region (B » Bcmb) yields the minimum 
gamma-ray flux, Fy ^in, in the hadronic model of radio mini- 
halos. We find Fy min to be a factor of 2 to 18 below the MAGIC 
ULs for spectral indices varying between 2.1 and 2.5. For a - 
2.2, following the universal CR model, the minimum gamma- 
ray flux is a factor of 3.2 lower than the MAGIC ULs. Second, 
by matching the radio emission profile (i.e., fixing the radial 
CR profile for a given magnetic field model) and by requiring 
the pion-decay gamma-ray flux to match the MAGIC flux ULs 
(i.e., fixing the normalization of the CR distribution) we obtain 
lower limits on the magnetic field distribution under considera- 
tion. Recall that we employ a parametrization of the magnetic 
field of B = Bo («/«{))""■ The inferred values range from 2 fiG < 
fio.min ^13 i-iG for the parameter space spanned by the magnetic 
field strength radial index, ob, and CR spectral index, a. Since 
fio.min is Smaller than recent field strengths estimates through 
Fara day RM studies in cool core clusters (e.g.. lKuchar & EnBlinI 
1201 1 ). this argues that the hadronic model is an interesting pos- 
sibility in explaining the radio mini-halo emission. This displays 
the potential of future gamma-ray observations of Perseus to fur- 
ther refine the parameters of the hadronic model and for eventu- 
ally assessing its validity in explaining radio (mini-)halos. This 
is true for the cuiTently operating Cherenkov telescopes and for 
the future planned Cherenkov Telescope AiTay whose sensitivity 
is planned to be about an order of magnitude higher than cuiTent 
instruments. 
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